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摘 要 运用 分 子 动力 学 方法 , 探究 bcc-W 在 中 子 辐 照 初期 , 由 辐 照 诱发 的 缺陷 形成 和 演化 的 微观 过 程 的 原子 机 制 . 选取 初 
AREJA T (PK A) BE Æ 1.0~25.0 keV, 模拟 温度 范围 100—900 K, 研究 W 中 位 移 级 联 产 生 的 缺陷 数量 及 分 布 , PKA 方 向 和 温 
度 对 稳定 Frenkel pairs 数 的 影响 , s [és AEE VA W Hf] s hz B fe. 结果 表明 , 车 级 联 诱 发 的 缺陷 在 峰值 阶段 呈 近 球形 密集 分 
布 , 稳定 阶段 Frenkel pairs 数 相 对 较 少 , 若 缺 陷 呈 非 球形 相对 分 散 分 布 , 稳定 阶段 Frenkel pairs 数 相对 较 多 ; 稳定 Frenkel 
pairs 数 受 PKA 方 向 的 影响 不 大 , 且 随 温度 升 高 有 下 降 趋势 ; 比较 而 言 , 间隙 原子 团 簇 分 数 比 空位 团 徐 分 数 高 , 而 空位 团 簇 倾 
向 于 形成 较 大 的 团 徐 ; W 的 平均 离 位 闵 能 受 温度 影响 较 小 , 并 具有 一 定 的 各 向 异性 . 
关键 词 W, 分 子 动力 学 , 位 移 级 联 , 稳定 Frenkel pairs Zt, BAKE 
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ABSTRACT Tungsten is a candidate material for the first wall and divertor in a tokamak fusion reactor, in 
which it is required to withstand a high neutron irradiation. The defects created in cascade form the primary state of 
damage and their subsequent evolution gives rise to important changes in their microstructures and engineering 
properties. However, the evolution and aggregation of radiation- induced defects in atomic level can not be ob- 
served by experiments up till now. In this work, molecular dynamics (MD) method was used to explore the micro- 
structural processes and atomic mechanism of the formation and evolution of defects in the initial stage of radi- 
ation in bcc- W. The range of primary knock-on atom (PKA) energies is 1.0—25.0 keV, and simulation temperature 
range from 100 to 900 K. The number and distribution of defects produced by displacement cascades have been 
studied; the influence of PKA direction and temperature on the number of steady Frenkel pairs has also been re- 
searched; defect clusters and the threshold energy have been simulated. The results showed that for morphology dis- 
tribution of defects induced in the peak time of cascade, the more intensive the defects are, the less the steady Fren- 
kel pairs numbers are, on the contrary, the more decentralized the defects are, the more the steady Frenkel pairs 


numbers are; the number of steady Frenkel pairs is insensitive to PKA direction, but has a trend to decline with 


* 


国家 


自然 科学 基金 项 目 21233010 和 51004012 资 助 


EE 


收 到 初稿 日 期 : 2014-09-04, 收 到 修改 稿 日 期 : 2015-03-27 
简介 : 姚 曼 , 女 , 1962 年 生 , 教授 , 博士 


DOI: 10.11900/0412.1961.2014.00492 


IL. 


LN 
IP 


202303.00613v1 


chinaXiv 


第 6 期 姚 曼 等 : W 辐 照 损伤 初 


ChinaXiv 合 作 期 刊 


期 的 分 子 动力 学 研究 725 


the temperature elevating; the percentage of interstitial clusters is higher than that of the vacancy clusters, while va- 


cancies tend to form larger clusters; the average threshold energy of W is less affected by temperature and has cer- 


tain anisotropy. The results of this work can provide data for analyzing the behavior of W material under nucle- 


ar environment. 


KEYWORDS W, molecular dynamics (MD), displacement cascade, steady Frenkel pairs number, displacement 


threshold energy 


无 论 在 International Thermonuclear Experimen- 
tal Reactor (ITER) Tokamak 或 未 来 的 Demonstration 
Fusion Reactor (DEMO) 聚 变 反 应 堆 中 , 第 一 壁 都 要 
承受 很 高 的 中 子 通 量 (2.5x102 n/em? RUP p AE TF 
高 的 热 通 量 (10 MW/m', ITER)". 因此 , 第 一 壁 的 保 
护 层 材料 , 所 谓 的 面向 等 离子 体 材料 , 必须 依据 此 材 


表面 的 辐 照 损伤 , 发 现 低 于 通常 认为 的 阀 能 也 能 引 
发 W 表 面 损伤 , 并 且 提 高 衬 底 温度 会 增加 W 原子 在 
表 


ri P^ AE XEBE RS n] Re TE, 这 种 表面 损伤 的 最 可 能 的 
能 量 范围 是 50.0~80.0 eV. Caturla 等 研究 了 W 在 
1.9 GeV 质子 辐 照 下 的 级 联 效率 和 缺陷 分 布 , 结果 
表明 , W 中 的 空位 是 独立 的 , 在 30.0 keV 的 PKA 作 


料 在 恶劣 辐射 条 件 下 的 性 能 进行 选择 . W 由 于 其 高 


热 导 率 、 低 溅 射 率 、 高 熔化 温度 ~ 和 低 毛 滞留 量 % 等 优 
点 成 为 未 来 聚变 堆 中 面向 等 离子 体 材料 的 领先 候 


| 下 没有 空位 团 簇生 成 . Troev 等 "通过 模拟 正 电子 
寿命 量子 机 制 计 算 来 研究 包含 和 He 的 W 中 的 缺 
陷 性 质 , HH 和 He 原子 被 吸收 在 W 的 空位 和 纳米 空 


选材 料 ; 偏 滤 器 可 以 清除 等 离子 体 中 的 HH 同 位 素 、 
He 和 杂质 , 与 第 一 壁 的 保护 层 一 样 , 它 也 倾向 于 由 
W 或 铝 合 金 制 成 ,并且 如 今 已 有 全 W 材 料 制作 的 偏 
滤器 ; 在 He 冷却 偏 滤 器 的 设计 中 , W 也 被 考虑 为 结 
构 材 料 外 . 因此 ,了解 W 在 辐 照 条 件 下 的 性 能 是 关乎 
核能 发 展 的 一 个 至 关 重 要 的 问题 . 但 是 , 辐 照 诱发 
点 缺陷 的 形成 与 聚集 发 生 在 皮 秒 (ps) 和 纳米 (nm) 时 
空 尺度 , 这 种 尺度 发 生 的 现象 目前 尚 没 有 实验 手段 
可 测 得 , 因此 , 计算 机 模拟 在 对 辐 照 效应 的 研究 中 
发 挥 了 关键 作用 外 其 中 , 分 子 动力 学 (molecular dy- 
namics, MD) 方 法 被 普遍 使 用 , 它 不 但 可 以 提供 原子 
动态 过 程 的 细节 , 并 可 以 将 结构 演化 作为 时 间 的 函 
数 , 而 且 , 它 能 描述 的 空间 和 时 间 尺 度 与 辐 照 引发 
的 位 移 级 联 的 空间 和 时 间 尺 度 一 致 . 

当 辐 照 中 子 进 入 材料 晶 格 当中 时 , 会 将 部 分 能 
量 转移 给 W 原子 , 此 W 原子 即 为 初始 碰撞 原子 (pri- 
mary knock-on atom, PKA), PKA 会 在 W 中 引发 一 
系列 的 碰撞 , 称 为 级 联 碰 撞 , 结果 产生 辐 照 位 移 损 
伤 , 形成 间 际 原子 和 空位 缺陷 对 (Frenkel pairs)、 空 位 
团 秘 、 间 际 原子 团 徐 以 及 位 错 环 、 空 位 四 面体 等 辐 


洞 中 , 而 了 再 和 He 在 较 大 的 纳米 空洞 中 的 存在 大 大 降 
低 了 正 电 子 寿 命 . Park 等 凸 研究 了 哑铃 型 间隙 原子 
的 分 布 方向 等 信息 , 得 出 稳定 阶段 的 哑铃 型 间隙 原 
子 都 沿 着 [111] 方 向 , 并 且 在 辐 照 损伤 过 程 中 , [111] 
方向 的 哑铃 型 间隙 原子 比 其 他 方向 更 快 与 空位 复 
合 . 随 着 W 材料 的 关注 度 不 断 提 升 , 人 们 需要 更 加 
深入 了 解 W 在 辐 照 条 件 下 的 缺陷 形成 机 制 和 各 种 
性 能 , 这 对 核 事 业 的 发 展 具有 重要 意义 . 

本 工作 利用 分 子 动力 学 方法 , 从 微观 原子 尺度 
上 系统 地 认识 W 在 中 子 辐 照 条 件 下 诱发 的 空位 、 间 
隙 原子 及 其 团 复 的 形成 机 制 , 计算 缺陷 数量 及 分 
布 , 分 析 温度 、 辐 照 能 量 `. PKA 入 射 方向 等 因素 的 影 
响 规律 . 另外 , 本 工作 还 计算 了 W B ES DL BU RE Ea. 
1 计算 方法 与 条 件 

W 有 cw 和 52 种 变型 . 在 标准 温度 0 "C (273.15 K) 
和 常 压 (1.01x10: Pa) F, a 型 W 是 稳定 的 bcc 结构. B 
型 W 只 有 在 有 和 氧 存 在 的 条 件 下 才能 出 现 , 它 在 
630 人 以 下 是 稳定 的 ,在 630 ‘CC 以 上 又 转化 为 a 型 
W, 并 且 这 一 过 程 是 不 可 逆 的 . 本 工作 研究 对 象 为 a 


照 缺 陷 . 在 辐 照 损伤 的 初期 阶段 , 级 联 碰撞 首先 会 


产生 大 量 的 Frenkel pairs, 并 在 短 时 间 内 达到 峰值 ， 
之 后 在 热 峰 的 作用 下 , Frenkel pairs 逐渐 减少 并 趋 于 


稳定 近年 来 , W 由 于 其 洪 在 的 优异 性 能 被 日 渐 重 
视 , Xu 等 "利用 MD 研究 了 辐 照 后 W 的 空位 扩散 行 
为 , REMEH W eM W 原子 反 冲 能 量 依赖 于 
PKA 方向, 并 且 反 冲 原 子 不 仪 与 最 近邻 空位 作用 ， 
也 会 与 第 二 近邻 的 空位 作用 . Li 等 利用 MD 研究 
不 同 温度 下 , 低能 量 (20.0~200.0 eV) He 原子 又 击 W 


型 W, Hl bec-W. 

本 工作 运用 分 子 动力 学 方法 研究 bcc-W 在 辐 照 
损伤 初期 缺陷 的 形成 和 演化 . PKA 能 量 Erca HW Js 
子 (质量 m2) 与 中 子 (能 量 Bo 质量 m1) 的 弹性 碰撞 公 
式 Es-AmmiPyEw (mi + mY f$ $1, 大致 关 系 为 Es 
46.05 Ew, 所 以 本 工作 研究 的 最 大 PKA 能 量 为 
25.0 keV, 大 约 对 应 1151.25 keV 中 子 能 量 . 聚变 堆 
中 14.0 MeV 中 子 能 量 大 约 对 应 306.0 keV 的 PKA 
AEE, 快 中 子 在 W 中 传递 给 PKA 原子 的 能 量 范围 
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在 1.0 到 100.0-200.0 keV 之 间 , 并 且 聚 变 材 料 产 生 
的 大 部 分 缺陷 都 是 小 于 100.0 keV 的 PKA 能量 的 作 
HE 根据 现 有 计算 条 件 , 本 工作 的 PKA 能 量 设 定 
73 1.0-25.0 keV, 而 针对 不 同 PKA 能 量 , 模拟 体系 
大 小 包含 约 20000~130000 个 原子 . 模拟 的 PKA 入 
射 方向 主要 为 [100] 和 [135]; 模拟 温度 为 100, 500 和 
900 K. 势 函 数 则 采用 ZBL 势 函 数 中 与 Lee 等 "提出 
的 2NN MEAM 势 函数 进行 而 合 , 从 而 能 够 更 好 地 
重 现 高 能 粒子 间 的 碰撞 问题 , 文献 [15] 中 也 采用 此 
类 型 势 函数 . 

首先 将 模型 在 设 定 好 目标 温度 的 等 温 等 压 系 综 
(NPT) F 381825 10 ps, 使 体系 能 量 达到 所 设 条 件 下 的 
最 低 状 态 , 获得 稳定 结构 (结构 1), 将 这 个 结构 作为 缺 
陷 分 析 的 参考 结构 . 然后 , 将 结构 1 中 某 一 原子 选 定 
为 PKA, 并 给 该 原子 按照 设 定 的 模拟 能 量 和 方向 重 
新 赋 一 个 初始 速度 , 其 他 原子 的 速度 和 方向 不 变 , 3 
将 此 结构 (结构 2) 作 为 辐 照 损伤 模拟 的 初始 结构 , 在 
正则 系 综 (NVT) 下 进行 一 定时 间 的 级 联 碰撞 模拟 . 一 
个 单独 级 联 的 模拟 时 间 取 决 于 模拟 体系 和 PKA 能 量 
的 大 小 以 及 势 函数 , 本 工作 模拟 发 生 在 10~30 ps 内 
的 原子 运动 行为 . 体系 外 部 采用 Berendsen 热 浴 来 将 
模拟 体系 的 温度 维持 在 本 工作 设 定 的 模拟 温度 , 整 
个 模拟 体系 采用 周期 性 边界 条 件 . 辐 照 结束 后 对 所 
得 结果 进行 统计 分 析 , 主要 通过 唱 格 匹配 分 析 方 法 
(lattice matching analysis, LMA)09 确 定 体 系 内 级 联 碰 
撞 诱发 点 缺陷 的 数目 及 分 布 情况 . 在 一 个 确定 的 能 
量 和 温度 的 条 件 下 进行 至 少 6 次 不 同 PKA 的 模拟 ， 
之 后 对 模拟 结果 进行 统计 平均 . 

离 位 阔 能 Es 即 转移 给 晶 格 原子 使 之 发 生 离 位 3 
产生 永久 缺陷 的 最 小 能 量 , 对 于 的 计算 , 本 工作 的 
模型 包含 4032 个 原子 , 共计 算 了 3 个 温度 (100, 500 和 
900 K), 3 个 方向 ([100], [110] 和 [111D 下 的 Es 在 确定 
一 个 温度 和 一 个 方向 的 条 件 下 , 对 选 定 的 PKA 原子 施 
加 一 定 区 间 的 能 量 , 此 区 间 大 小 为 15 eV, 以 1 eV 3$ 
增 , 在 此 能 量 区 间 , 将 产生 一 个 稳定 Frenkel pair 所 
对 应 的 最 小 PKA 能 量 定义 为 该 原子 在 所 选 方 向 上 
的 ,这 种 模拟 共 进行 75 次 , 每 次 选取 不 同 的 PKA， 
之 后 将 75 次 计算 取 平 均值 作为 W 在 该 温度 和 该 方 
向 下 的 Es. 其 他 温度 和 方向 的 计算 也 如 上 所 述 , 对 
的 计算 总 次 数 为 675 次 . 

2 计算 结果 与 分 析 
2.1 W 辐 照 过 程 的 Frenkel pairs 随时 间 演 化 及 形态 
分 布 
1 所 示 为 W 的 Frenkel pairs 在 不 同 PKA 能 量 
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—e— 5.0 keV 
—e— 15.0 keV 
—A— 25.0 keV 


Number of Frenkel pair 
n 
N 
[«] 
e 


1 100 开 时 不 同 初始 碰撞 原子 (PKA) 能 量 下 辐 照 诱发 
Frenkel pairs 数量 随时 间 的 演化 曲线 


Fig.1 Evolution curves of number of Frenkel pairs with 


time in displacement cascade casued by a set of pri- 
mary knock-on atoms (PKAs) with various energies 
at 100 K 


(5.0, 15.0 4125.0 keV) 下 随时 间 的 演化 曲线 . 曲线 的 
特征 与 之 前 研究 结果 3 一致. 在 初始 时 刻 , Fren- 
kel pairs 迅速 增加 并 达到 峰值 (弹道 相 阶段 ), 然后 快 
速 下 降 ( 回 复 阶 段 ), 之 后 到 达 稳 定 阶段 (10~30 ps). 
可 以 看 出 , 达到 峰值 和 稳定 状态 的 时 间 都 随 PKA 能 
量 (1.0~25.0 keV) 的 增加 而 增加 . 
图 2 所 示 为 PKA 能 量 为 10.0 keV 时 , 点 缺陷 分 
别 在 峰值 时 刻 和 稳定 阶段 的 Frenkel pairs 分 布 图 . 
黑色 原子 代表 空位 , 主要 分 布 在 级 联 中 心 处 ; 白色 
原子 代表 间 隐 原子 , 多数 分 布 在 级 联 中 心 的 边缘 . 
可 以 看 出 , 级 联 碰撞 经 过 峰值 阶段 之 后 , 在 热 峰 的 
作用 下 大 量 间 辽 原子 与 空位 发 生 复合 重组 , 使 Fren- 
kel pairs 数 迅 速 下 降 , 最 终 在 稳定 阶段 留 下 较 少 的 
Frenkel pairs. 

模拟 结果 显示 , 从 能 量 超过 10.0 keV 开始 , 不 同 
PKA 引 发 的 缺陷 峰值 阶段 的 级 联 分 布 形态 变化 程 
度 变 大 , 这 会 影响 最 终 Frenkel pairs 的 数量 和 分 布 ， 
如 图 3 所 示 . 总 体 来 讲 , 峰值 缺陷 形态 大 致 分 为 2 
种 : 近 球 形 密集 型 和 非 球形 相对 分 散 型 . 如 果 峰 值 
时 刻 缺陷 呈现 近似 球形 的 密集 分 布 , 那么 通常 峰值 
缺陷 数目 比较 多 , 由 于 空位 和 间 隐 原子 的 分 离 距离 
较 近 , 热 峰 作用 时 间 较 长 , 使 缺陷 复合 相对 容易 , 导 
致 稳定 阶段 Frenkel pairs 相对 较 少 的 几率 大 , 如 图 
3a 所 示 ; 若 峰 值 时 刻 缺陷 呈现 非 球形 的 相对 分 散 分 
布 , 那么 峰值 缺陷 数目 相对 少 , 热 峰 阶段 也 会 相对 
短 , 且 间 隙 原子 与 空位 分 离 距离 较 远 , 不 利于 缺陷 
的 复合 , 导致 残留 较 多 的 Frenkel pairs, 如 图 3b 所 
示 ; 且 这 种 趋势 受 PKA 方 向 的 影响 不 明显 . 

峰值 时 刻 缺 陷 为 非 球形 相对 分 散 型 时 , 有 时 还 
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2100 K 下 PKA 能量 为 10.0 keV, PKA 沿 [135] 唱 向 的 级 联 碰撞 中 不 同时 刻 Frenkel pairs 的 分 布 情况 
Fig.2 Distributions of radiation-induced Frenkel pairs in different moments at 100 K (The PKA energy is 10.0 keV and the 


PKA velocity direction is parallel to the crystallographic direction [135] of bcc- W, black balls represent the vacan- 
cies and white balls the interstitials) 


(a) peak time 1.04 ps (b) steady stage 12.22 ps 


ui 
四 | 
EPE 
(a) 


3 100 KK 下 PKA 能 量 为 20.0 keV 时 峰值 时 刻 不 同 形态 下 的 Frenkel pairs 的 分 布 情况 
Fig.3 Distributions of radiation-induced Frenkel pairs at peak time in different morphologies at 100 K (The PKA energy is 


(b) 


20.0 keV, the PKA velocity direction is parallel to the crystallographic direction [135] of bcc-W, and in the circles 
the alternating sequence of interstitials and vacancies are shown) 
(a) dense morphologies, the number of Frenkel pairs in peak time: 2631, the number of Frenkel pairs in steady stage: 
23 
(b) dispersed morphologies, the number of Frenkel pairs in peak time: 829, the number of Frenkel pairs in steady 


stage: 32 


可 观察 到 沿 密 排 [111] 方 向 的 空位 和 间隙 原子 交 蔡 文献 [3] (a=3.706, b=0.590) 和 文献 [9] (a=3.335, b= 
序列 , 如 图 3b 中 圆圈 部 分 所 示 , 间隙 原子 通过 这 种 ”0.602) 等 的 研究 结果 相似 . 图 4 也 显示 了 利用 Nor- 
形式 迁移 到 远离 级 联 中 心 的 位 置 , 使 此 部 分 间隙 原 ” gett-Robinson-Torrens (NRT) 公 式 "" 来 估算 W 辐 照 后 


子 不 易 被 复合 , 也 会 导致 Frenkel pairs ££ €. 的 Frenkel pairs 个 数 Nrswrn: 
2.2 稳定 阶段 Frenkel pairs 数 人 = 0.8E,, / QE) (2) 


本 工作 计算 的 稳定 阶段 Frenkel pairs 数 Ns 与 。 式 中 , HUE Enos, 尽管 Eww 是 指 损伤 能 量 而 不 是 
PKA 能 量 的 关系 如 图 4 所 示 , 符合 Bacon 等 "推荐 。 PKA 能 量 , 但 MD 模拟 不 考虑 电离 能 损 , 所 以 用 Eva 


= 


的 经 验 公式 : 近似 代 闪 ,来 比较 MD 模拟 所 得 的 Ns 预测 什 是 合 
Ny = aE’ka (1) x&pgusm. 选取 为 85 eV ^91, 
式 中 ,a 和 是 受 辐 照 材料 的 种 类 和 温度 影响 较 小 的 由 图 4 可 知 , 与 NRT 预测 值 相 比 , MD 模拟 的 级 


常量 , 将 MD 数据 进行 拟 合 后 , a=3.433, b=0.697, 与 联 产 生 的 Nes 显 著 降 低 , 这 是 因为 级 联 产 生 的 间 际 原 
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Fig.4 Numbers of survival Frenkel pairs at steady stage as 
a function of PKA energy (Solid dots represent the 
numbers of steady Frenkel pairs estimated by Nor- 
gett- Robinson- Torrens (NRT) formula Nr», the 
hollow circles are the values simulated by molecular 
dynamics (MD) Nrs, the bars represent the standard 


deviations) 


子 只 有 很 小 一 部 分 是 通过 聚焦 碰撞 链 的 形式 迁移 
到 远 处 , 而 大 多 数 的 间隙 原子 产生 于 级 联 不 规则 区 
域 的 边缘 , 这 些 间 际 原 子 与 缺陷 中 心 的 空位 距离 非 
常 近 , 并 且 热 峰 阶段 级 联 缺 陷 中 心 有 很 高 的 动能 ， 
有 助 于 间 际 原子 和 空位 的 复合 , 以 致 产生 的 Frenkel 
pairs 显著 减少 , 文献 [12] 计 算 的 W 和 文献 [17] 中 介 
绍 的 Fe, Cu, Zr 等 多 种 金属 得 到 相 类 似 的 结果 . 由 这 
2 种 估计 材料 辐 照 损伤 缺陷 的 方法 引出 了 级 联 效率 . 

级 联 效 率 定 义 为 mFNisy/Nrsaro". 随 着 PKA 能 
量 的 增 大 , 级 联 效率 呈 降 低 趋势 , 最 终 趋 于 稳定 ， 
约 为 0.3 (图 5), 与 文献 [19] 所 得 值 相符 , Bacon 等 四 
的 研究 也 表明 , 不 同 金属 的 级 联 效率 通常 稳定 在 
0.2~0.4. 
2.3 Nes 与 PKA 入 射 方向 关系 

众所周知 , 晶体 中 不 同 唱 向 原子 的 排列 周期 及 
ij 95 Fe BE — XE 28 ke, 致使 其 性 能 具有 各 向 异性 ， 
因此 , 在 级 联 碰撞 中 bcc-W 也 将 显示 出 辐 照 性 能 的 
各 向 异性 . 表 1 为 W 在 不 同方 向 的 PKA 作用 下 , 稳 
定 阶 段 体 系 中 的 Frenkel pairs 数 和 Frenkel pairs 数 
随 方向 的 浮动 情况 . 

Frenkel pairs 个 数 的 浮动 采用 以 下 公式 计算 : 

F - 10096 x (N, -NY N.. Q3) 


IP, F V Frenkel pairs 数 随 不 同 PKA 方 向 的 浮 
动 , Y 代 表 每 个 方向 的 稳定 Frenkel pairs 数 , Ne 代表 
5 个 方向 Frenkel pairs 数 的 平均 值 . 结果 表明 , 数据 浮 
动 最 大 只 有 10% 左 右 , 如 表 1 所 示 . 可 以 看 出 , 本 工作 
研究 的 5 个 方向 的 稳定 Frenkel pairs 个 数 在 3.0 keV 
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0.2 
0 
PKA energy / keV 
5 级 联 效 率 n 随 PKA 能 量 的 变化 
Fig.5 Cascade efficiency 7 as a function of PKA energy 


表 1 100 K, Ev 733.0 keV 时 不 同 PKA 方 向 下 单 晶 W 中 

级 联 碰撞 诱发 稳定 Frenkel pairs 数 和 浮动 情况 
Tablel Number and its floating of Frenkel pairs at steady 
stage along various crystallographic directions of 


bcc-W with PKA energy of 3.0 keV at 100 K 


PKA direction Nis F 
[100] 7.00+0.02 —4.4% 
[110] 6.90+0.22 -5.8% 
[111] 6.67+0.02 -8.9% 
[112] 7.97+0.04 8.9% 
[135] 8.07+0.14 10.2% 
~ Note: M«—number of Frenkel pairs, F—floating of — 
Frenkel pairs 
能 量 下 相差 很 小 


所 以 在 低能 量 下 , PKA 方 向 对 bcc-W 级 联 碰撞 
稳定 阶段 Frenkel pairs 数 没 有 明显 影响 . 文献 [21] 
在 PKA 能 量 为 0.3~5.0 keV 对 hcp-Ti 进 行 的 计算 研 
2t, 以 及 文献 [16] 均 获得 类 似 结 果 . 需要 指出 的 是 ， 
为 了 增强 模拟 结果 的 准确 性 , 本 工作 在 每 个 方向 下 
进行 至 少 3 次 不 同 PKA 的 模拟 , 对 模拟 结果 进行 统 
计 平 均 . 

2.4 Nis 与 温度 的 关系 

为 研究 温度 对 辐 照 缺陷 产生 的 影响 , 本 工作 在 
PKA 能 量 3.0 keV, 温度 分 别 为 100, 500 #1900 K F 
进行 模拟 , 每 个 温度 下 进行 9 次 计算 , 结果 如 图 6a 所 
示 . 可 以 看 出 , Nis 随 温度 的 升 高 有 轻微 下 降 . 为 了 在 
不 同 能 量 下 验证 , 又 计算 了 5.0, 10.0 以 及 15.0 keV 
下 的 情况 , 如 图 6b 所 示 . 总 体 来 看 , Nis 随 着 温度 的 
升 高 大 致 呈 略 微 下 降 趋势 , 与 文献 [12,22,23] 得 到 的 
结果 相似 , 下 降幅 度 在 7%~40% 之 间 . 这 可 能 是 因为 
温度 升 高 的 同时 增加 了 W 原子 活动 性 与 热 峰 持续 
时 间 , 两 者 在 一 定 程 度 上 促进 了 间隙 原子 与 空位 的 
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6 Nes 随 温度 的 变化 
Fig.6 N« at steady stage in cascade simulations with PKA energies of 3.0 keV (a) and 5.0~15.0 keV (b) at various tempera- 


tures 100-900 K of bcc- W 


复合 , 使 更 多 缺陷 复合 漂 灭 , 从 而 减少 稳定 阶段 留 
存 的 Frenkel pairs Zi. 男 外 一 个 可 能 的 原因 是 温度 
升 高 使 级 联 中 的 换 位 碰撞 序列 (replacement colli- 
sion sequence, RCS) 长 度 变 短 , 导致 级 联 趋向 于 以 紧 
密 的 形式 存在 , 从 而 使 空位 和 间隙 原子 分 离 的 距离 
较 小 , 有 助 于 复合 作用 的 发 生 . 
2.5 缺陷 团 簇 
缺陷 团 簇 影响 点 缺陷 的 迁移 方向 和 难 易 , 还 可 
影响 点 缺陷 在 级 联 中 的 回复 和 洒 灭 以 及 材料 辐 
照 后 的 性 能 演化 , 所 以 团 艇 的 作用 不 可 小 视 . BLA 
} 数 矿 是 指 形成 缺陷 团 秘 的 缺陷 个 数 入 与 此 类 型 
缺陷 的 总 的 缺陷 数 Nu 的 比值 : 
六 EM。 (4) 
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Fig.7 Fractions of interstitials fa; and vacancies fa, via PKA 


energy (10.0—25.0 keV) at 100 K at steady stage 


例如 , EBR JE T BIE B fes JFE Yo. 8] B Je. ^ H1] 
TE BI PEI BR JE ^ I 2. Nu 与 间隙 原子 总 个 数 NMu 的 比 
值 . 图 7 所 示 为 100 K 下 , 稳定 阶段 记 与 空位 团 簇 分 
Ji fa BE PK A BETIS AE 
H E] 7 n] AD, 缺陷 团 簇 分 数 大 致 在 0.15~0.40 之 
间 , 与 文献 [12,19,22] 相 似 , 团 簇 分 数 不 高 ,可 见 , W 
中 大 部 分 缺陷 仍 以 单 缺陷 的 形式 存在 , 但 相 比 之 
下 , 大 部 分 情况 正高 于 fs 与 文献 [3,12,19] 结 果 一 
致 . 以 上 研究 表明 , 单个 间隙 原子 或 间隙 原子 团 簇 中 
大 部 分 原子 在 MD 模拟 的 时 间 尺 度 内 都 具有 很 高 的 活 
动 性 . 据 文献 报道 W 原子 的 迁移 激活 能 在 [111] 方 向 为 
(0.03 + 0.002) eV**, (0.047 2 0.01) eVP, 0.03 eV 忠和 
0.061 eV??, 而 W 的 空位 迁移 能 为 1.44 eV, 1.78 eVP9 
和 2.054 eVea, 室 位 平均 迁移 能 为 (3.1+0.05) eV”, 
由 此 可 见 , 间隙 原子 迁移 能 要 远 远 小 于 空位 迁移 
能 , 而 迁移 能 小 使 迁移 更 加 容易 , 从 而 更 易 聚 集 形 
RAI, 并 且 已 有 研究 表明 , 空位 一 旦 形成 便 不 易 
移动 , 例如 Fikar 等 中 在 1000 K 时 , 在 一 个 大 体系 中 


jami 


引进 1 Ô Frenkel pair, MD 模拟 进行 50 ps 时, MJ [f 
原子 沿 [111] 方 向 迁移 15 个 晶 格 距离 , 空位 没有 发 
生 移动 ; Bai $E 4 Cu 的 研究 也 得 出 空位 不 易 移 动 
的 结论 . 本 工作 也 对 级 联 过 程 稳 定 阶段 进行 了 可 视 
化 , 发 现 W 的 空位 形成 后 基本 不 发 生 移动 (图 8). 所 
以 间隙 原子 可 能 会 凭借 迁移 较 易 的 优势 更 易 形成 
zip 

本 工作 计算 结果 还 显示 , 在 低能 量 (小 于 5.0 keV) 
PKA 的 作用 下 , 无 论 是 间隙 原子 还 是 空位 都 很 难 形 
成 团 徐 . PKA 能 量 增 大 时 , 单个 缺陷 仍 占 总 缺陷 数 
的 很 大 一 部 分 比例 , 而 由 图 9 可知 , 与 间隙 原子 团 簇 
相 比 , 虽然 空位 团 徐 大 小 ( 即 一 个 空位 团 簇 中 包含 空 
位 的 个 数 ) 为 4 个 或 小 于 4 个 的 相对 少 , 但 其 更 易 形 
成 大 团 簇 ,形成 空位 个 数 大 于 5 的 团 簇 , 甚至 形成 一 个 空 
位 个 数 为 12 的 大 团 复 有 趣 的 是 , 对 于 W 的 空位 团 簇 的 
研究 结果 有 很 大 分 歧 , 文献 [3] 在 PKA 能 量 为 20.0 keV, 
级 联 进行 5 ps 时 出 现 空位 团 簇 ,但 在 30 ps 时 团 簇 又 消 
失 了 ; 质子 辐 照 W 的 研究 结果 "表明, W 的 空位 是 独 
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8 900 K 下 PKA 能 量 为 15.0 keV 时 稳定 阶段 点 缺陷 (间隙 原子 和 空位 ) 随 时 间 演 化 
Fig.8 Evolution of point defects during steady stage at 900 K with PKA energy 15.0 keV at 20.52 ps (a), 25.90 ps (b), 
28.83 ps (c), 31.76 ps (d), 34.70 ps (e) and 40.72 ps (f) (Black balls represent the vacancies and white balls the inter- 
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Fig.9 Relationship of cluster size (the number of defects in 


each cluster) and number with PKA energy 25.0 keV 
at 100K 


立 的 , PKA 能 量 为 30.0 keV IRE UTE 254 AEE pic; 
日 文献 [22] 研 究 得 出 , 10.0 keV 的 PKA 能 量 就 能 
致 大 的 空位 团 簇生 成 . 所 以 , 关于 W 辐 照 下 空位 团 
簇 问题 还 需 深入 研究 . 
2.6 BBE E, 

实验 和 模拟 都 对 W AI Est T RUE BRL E 76, 但 结 
果 差 别 很 大 . 并 且 实 验 中 只 可 以 测量 到 缺陷 形成 的 
继 发 效应 (电阻 或 位 错 环 ), 因此 真实 很 难 通过 实 
验方 法 测 得 . 分 子 动 力学 模拟 是 一 种 计算 的 理想 


= 


方法 , 它 可 以 直接 “控制 ?特定 原子 在 唱 格 特定 方向 
上 的 单个 位 移 事 件 , 不 受 其 他 事件 的 影响 . 
图 10 为 在 [100], [110] 和 [111] 3 个 方向 , 温度 为 
100, 500 和 900 K 的 条 件 下 对 W B Ei S £s TR, B 
中 还 引用 其 他 学 者 的 研究 工作 ”2 作为 比较 . 可 
以 看 出 , 本 工作 模拟 结果 基本 在 先前 研究 范围 2 
内 , 在 100, 500 和 900 下 的 的 平均 值 (3 个 方向 
取 平 均值 ) 分 别 为 62.17, 62.77 和 59.47 eV. 可 见 , Ea 
随 温度 的 变化 幅度 不 大 , 到 900 K 时 有 了 上 略微 下 降 ; 但 
从 单个 方向 上 看 , 密 排 [111] 方 向 的 包 随 温度 升 高 而 
略微 上 升 , [100] 和 [110] 方 向 的 五 略微 降低 . 

随 温度 升 高 而 上 升 可 能 是 由 于 以 下 原因 : 体 
系 中 缺陷 的 产生 主要 是 由 于 级 联 碰撞 链 以 及 碰撞 
后 容易 在 密 排 方向 产生 的 RCS. 高 温 加 剧 了 体系 中 
原子 的 热 振 动 , 使 RCS 长 度 变 短 , 缩短 了 空位 和 间 
隙 原子 的 平均 分 离 距离 , 使 复合 的 可 能 性 更 高 , 同 
时 除 RCS 以 外 的 间 际 原子 与 空位 的 回复 作用 也 得 
到 促进 , 所 以 需要 更 多 的 能 量 去 产生 稳定 的 Frenkel 
pairs, 从 而 使 轧 升 高 . Bs 随 温 度 升 高 而 降低 则 可 能 是 
于 高 温 不 仅 促进 了 回复 , 同时 也 给 某 些 间隙 原子 
提供 了 通过 热 激 活 迁 移 到 远离 级 联 中 心 的 位 置 的 
机 会 , 从 而 导致 的 降低 . 

基于 以 上 分 析 , 温度 对 的 影响 是 以 上 2 种 效 
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果 相 反 的 机 制 共同 作用 的 结果 , 哪个 机 制 占有 更 大 
的 优势 , 5 就 会 表现 出 相应 的 温度 效应 . 

图 11 为 与 PKA 入 射 方向 关系 的 盒 形 图 . 可 见 ， 
PKA 入 射 方向 对 E 的 影响 很 大 , 这 与 之 前 许多 研究 结 
果 一 致 . ,在 [110] 方 向 平均 值 呈 现 最 大 值 (80.0 eV), 
在 密 排 方向 [111] 平 均值 呈现 最 小 值 (40.0 eV), 密 排 
方向 相对 小 可 能 是 由 于 沿 着 此 方向 的 RCS 使 间 
隙 原子 更 容易 迁移 到 离 空 位 较 远 的 地 方 , 产生 更 多 
Frenkel pairs. 文献 [30] 对 Cu 和 文献 [31] 对 V 的 计算 
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Vi, 在 相同 的 辐 照 能 量 ( 几 千 电 子 伏 ) 下 , 低 ELI AE TE 
会 使 体系 中 产生 更 多 缺陷 . 这 与 前 文 所 述 的 稳定 
Frenkel pairs 数 随 温度 升 高 呈 降 低 趋势 似乎 相 了 矛盾 ; 
H. Ei 受 一 定 的 方向 影响 , 这 与 前 文 所 述 的 稳定 Fren- 
kel pairs Zi fE PKA 能 量 为 3.0 keV 时 对 PKA 方 向 的 
敏感 度 小 相 了 矛盾 . 这 是 因为 , 级 联 碰撞 过 程 除 了 原 
子 之 间 的 碰撞 移 位 以 外 , 相伴 的 热 峰 效应 不 可 忽 
各 , 热 峰 会 对 点 缺陷 的 回复 产生 很 大 影响 , 进而 影 
响 辐 照 诱发 的 稳定 Frenkel pairs, MD 模拟 的 级 联 碰 


结果 也 表明 , 密 排 方向 的 相对 其 他 方向 较 低 ; 本 
工作 计算 的 [100], [110] 和 [111] 方 向 均 为 低 指数 方 
I], JC I] EZ] 73 90.0 eV, 文献 [32] 对 Cu 的 计算 结 
TE IK 8 207] E ES SEMIS, [100] 和 [110] 方 向 只 有 
25.0 eV 左右 , 而 高 指数 方向 及 相对 较 高 , 可 达到 
180.0 eV. 可 见 , 具有 一 定 的 各 向 异性 . 

可 见 , 温度 (100~900 KK) 对 二 平均 值 的 影响 不 太 
明显 , 到 900 K ET, 严 略 微 降低 3.0 eV 左右 , Tfi E DE 
低 说 明 产 生 稳 定 缺 陷 所 需要 的 能 量变 小 , 也 就 是 
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Fig.10 Comparisons of threshold energies simulated in this 


work and other works via simulations or experi- 


ments for W 


100 


J 


[100] [110] [111] 
PKA direction 


Threshold energy / eV 


Color online 
图 11 100 K FARI PKA Zr E FIEL Ea 
Fig.11 Threshold energy at 100 K along different crystallo- 
graphic directions (Blue dots represent average val- 
ue and red dots represent the values with the fre- 


quency of maximum) 


撞 产 生 的 Frenkel pairs 数 远 小 于 NRT 预测 值 , 正 是 
热 峰 作用 的 一 个 典型 体现 , 而 在 离 位 阔 能 的 计算 当 
中 , 由 于 施加 的 PKA 能 量 很 小 ( 几 十 电子 伏 ), 所 以 热 
峰 的 作用 很 小 , 以 至 于 温度 对 的 综合 作用 主要 体 
现在 使 缺陷 对 分 离 , 从 而 导致 b 随 温度 升 高 呈 略 微 
下 降 趋势 . 方向 对 的 作用 则 因为 晶 格 排列 的 周期 
及 踢 密 不 同 而 呈现 出 一 定 的 各 向 异性 . 因此 , 尽管 
是 影响 间 孙 和 空位 缺陷 形成 的 重要 参数 , 但 是 最 
终 的 稳定 Frenkel pairs 数 不 能 简单 的 单 靠 的 大 
小 来 判断 . 

3 结论 

(1) 在 对 bcc-W 模 拟 的 PKA 能 量 (1.0~25.0 keV), 
模拟 温度 100—900 K F, 级 联 峰 值 阶段 缺陷 分 布 大 
致 呈现 2 种 形态 : 近 球 形 密集 型 和 非 球形 相对 分 散 
型 . 多 数 情况 下 , 峰值 缺陷 呈 近 球形 密集 型 分 布 时， 
稳定 Frenkel pairs 数 相对 较 少 , 反之 呈 非 球形 相对 
分 散 型 分 布 时 , Frenkel pairs 数 相 对 较 多 . 

(2) 稳定 Frenkel pairs 数 在 PKA 能 量 为 3.0 keV, 
温度 为 100 K 时 受 PKA 入 射 方向 (本 工作 计算 
[100], [110], [111], [112] 和 [135]) 的 影响 不 大 , 各 方 
向 稳定 Frenkel pairs 数 的 浮动 最 大 在 10% 左 右 ; 
PKA 能 量 为 3.0~15.0 keV F, 稳定 Frenkel pairs 数 随 
温度 (100~900 KR) 升 高 而 略微 降低 . 

(3) Œ PKA (1.0-25.0 keV) 作 用 下 , W 中 缺陷 多 
以 单 缺 陷 形式 存在 , 难于 形成 缺陷 团 复 , HER) 
多 数 在 0.15~0.40 之 间 , 比较 而 言 , 间隙 原子 团 簇 分 
数 高 于 空位 团 艇 分数 ; 在 PKA 为 25.0 keV 时 , 单个 
团 徐 中 缺陷 数 为 4 个 或 少 于 4 的 间隙 原子 团 徐 要 多 
于 空位 团 簇 , 而 单个 团 簇 中 缺陷 数 多 于 5 个 的 团 簇 
只 有 空位 团 簇 , 空位 更 倾向 于 形成 大 团 簇 . 

(4) 离 位 阔 能 受 温度 影响 较 小 , 在 100, 500 和 
900 KK 下 3 个 方向 [100], [110] 和 [111] 的 离 位 闵 能 都 在 
60.0 eV 左右 ; 离 位 闵 能 呈现 一 定 的 各 向 异性 , 100 K 
时 在 密 排 方向 [111] 的 最 小 , 平均 值 为 40.0 eV, 在 
[110] 方 向 的 最 大 , 平均 值 为 80.0 eV. 
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